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Abstract 

In the present paper we consider general counting processes stopped at a 
random time T, independent of the process. Provided that T has the decreas- 
ing failure rate (DFR) property, we give sufficient conditions on the arrival 
times so that the number of events occurring before T preserves the DFR 
property of T. These conditions involve the study of the conditional interar- 
rival times. As a main application, we prove the DFR property in a context 
of maintenance models in reliability, by the consideration of Kijima type I 
virtual age models under quite general assumptions. 
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1 Introduction 



Decreasing failure rate (DFR) is a property describing a system which improves with 
age. For instance, the lifetime of a mechanism whose distribution is a mixture of 
exponential random variables has this property (see Proschan |20j). In fact, it is 
well known that mixtures of DFR distribution have the DFR property (see also [H] 
or [T5], for instance, for different contexts in which the DFR property can arise). 
Motivated by the fact that the DFR property is easily preserved by mixtures, and 
with the aim of completing the previous results obtained in j6], we will study the 
discrete DFR property for N(T), in which {N(t) : t > 0} is a counting process and 
T is a random time independent of this process. The discrete DFR property is a 
logconvexity condition for the survival function (see Esary et al. [10J or Grandell 
H31 Ch. 7]), that is, 

P{N(T) > n + l) 2 < P(N(T) > n)P(N(T) > n + 2) n = 0, 1, . . . . (1) 

Note that in this case we consider a mixture of random variables which, in 
general don't have the DFR property (for instance, in a Poisson process, each N(t) 
is discrete increasing failure rate). 

The interest in the study of N(T) comes back to |9J, p. 42 in which several ex- 
amples were given. Natural applications arise also in queuing theory, when studying 
the stationary number of waiting customers (see [2], for instance). Recent applica- 
tions have been found in health sciences (see [22] ), in which the counting process 
describes the number of tests for a disease (HIVS for instance) in a person at risk 
until the random time T at which this person is infected. 

Preservation of different ageing properties for a classical renewal process (that 
is, with independent, and identically distributed interarrival times) have been ex- 
tensively studied (see JTUJ, [T31 EJ 0], for instance). Recently, attention is focused in 
generalizing these properties for counting processes exhibiting dependence between 
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the interarrival periods (see [2TJ |3]). These models are of interest in reliability, for 
instance, when considering a process which models the successive repairments of a 
mechanism which deteriorates with age (see p2]). In the study of these models, 
it is sometimes usual to consider the simplifying hypothesis that the interarrival 
epochs are discrete or absolutely continuous random vectors and one of the aims of 
this paper is to present our results in a general background easily adapted to usual 
examples. 

As mentioned before, we complete the DFR preservation results obtained in 
[6] , going deeper in the study of the dependence properties between the interarrival 
times which can give place to this property, with special attention to get preservation 
results for general counting models of interest in reliability. Our main result in this 
paper is Proposition I3.4[ which shows the preservation of the DFR property under 
a (strong) decreasing behaviour of the interarrival epochs. The result is based on 
Propositions 13.11 and 13.21 The last one provides a way to check (CQ) for n = 1, 
when the second renewal epoch satisfy our decreasing condition. Proposition 13.11 
provides a way to check ([1]) for a general n by means of the conditional distribution 
of (X n+ i, X n+2 ) given the past. Our result differs from the ones given in [6], in 
the sense that the last ones need, along with (different) decreasing conditions for 
the renewal epochs, some additional conditions. As an application we will show 
in Section 3 the DFR property for Kijima type I models (see [17]) under more 
general assumptions than in 0, Section 5]. Finally, as an immediate application of 
Proposition 13 .21 we will give another preservation result which assumes the condition 
used in [6J Proposition 3.1. (a)]) for the first interarrival epochs, together with the 
same condition for the the next two renewal epochs knowing the past history. 

The proofs of our results are made without considering the absolutely continuity 
assumption on the interarrival epochs before mentioned, so that in Section [2] we 
introduce some preliminar background in order to specify the conditional distribu- 
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tion for two inter renewal times given the past history, making use of conditional 
probability kernels. Also in Section [2j we give some preliminary definitions which 
are going to be used along the paper. 

2 Preliminaries 

Let {N(t) : t > 0} be a general counting process whose arrival times are denoted by 
S n , n = 0, 1, . . . (5*o = 0) and the interarrival epochs are denoted by X n , n = 1, 2, . . ., 
that is: 

n 

S n = ^2Xi, n=l,2,... 

i=l 

The only assumption about the Xi is that they are non negative random variables. 
Recall that the counting process is defined through the renewal epochs by the fol- 
lowing expression 

N(t) = max{n : S n <t}, t > 0. 

Our paper deals with the preservation of the decreasing failure rate property. 
This concept, for discrete random variables, was reminded in the Introduction. The 
concept for 'general' random variables is recalled below (see Barlow and Proschan 
[7j) along with another one reliability class (New Worse than Used), which is weaker 
than the DFR property, and will be used throughout the paper. 

Definition 2.1. Let X be a nonnegative random variable with G and G := 1 — G 
the corresponding distribution and reliability functions. X (or G) is said to be: 

1. Decreasing failure rate if 

G(z + t)/G{t) is increasing in t, for all z > 0; 

2. New worse than used (NWU) if 
G(z)G(t) < G{z + t), for all z, t > 0; 
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Note that if the monotonicity and sense of the inequalities in Definition 12.11 are 
reversed, we obtain the dual concepts of increasing failure rate (IFR) and new better 
than used. 

Our results make use of comparisons of random variables in the usual stochastic 
order, which we recall in the following. 

Definition 2.2. Let X and Y be two random variables with distribution functions 
Fx and Fy , respectively. X is said to be smaller than Y in the usual stochastic order 
(denoted X <$t Y) if Fx it) > Fy(t) for all realt or, equivalently, if Fx(t) < Fy(t) 
for all real t. 

In our results concerning the DFR preservation property in a renewal process, 
we will make use of the conditional distributions of the inter renewal epochs. More 
specifically, for a given natural number m, we will consider the conditional distri- 
bution of (X n+ i, . . . , X n+m ) given that (X 1; . . . X n ) is known. In order to deal in a 
general setting, we will use conditional probability kernels. Consider two Borel sets 
SCR" and T C R m . Denote by B(T) the family of Borel sets on T. Recall that a 
probability kernel [i from S to T is a mapping 

At : S x B(T) -> [0, 1] 
(x,£) ~» ^(B) 

verifying the following properties: 

1. /i'(.B) is a measurable function for B e B(T) fixed. 

2. yU x (- ) is a probability distribution on (T, B(T)), for x G S fixed. 

Let (Xi, . . . , X n , X n+1 , . . . X n+m ) be a nonnegative random vector defined on a given 
probability space. Taking into account [121 Thm. 6.3, p. 107]), the existence of 
a probability kernel [i from R™ to R™ such that for each Borel set B C R™, 
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l^ Xl '-' Xn \B) is a version of P((X n+ i, . . . ,X n+m ) G B\(Xi, . . . ,X n )) can be always 
guaranteed. Moreover, if we denote by F n the distribution function of (Xi, . . . , X n ), 
we have for a nonnegative measurable function / : R™ — > R + that (cf. [121 p-108]) 

E[f(X u X n+m )} = [ /(z)rfF n+m (z) = / dF n (x) / /(x, y)^ x (y) 

-/zGR™ +m ./xGR™ </ yGR™ 

and obviously, if we take a set iV n G i3(R") such that P((X 1; . . . ,X n ) G iV n ) = 1, 
we have 

E[f(X 1 ,...,X n+m )}= [ dF n (x) / /(x,y)^ x (y). (2) 

ixeJV„ ./ye Rip 

In fact, (J2J) is the property we will need in our proofs. Observe that (j2J) is satisfied 
for all nonnegative / if an only if it is satisfied for /(x, y) = 1axb( x > y)> f° r Borel 
sets A C jV n , 5 C R™, in which 1ax_b(' ) denotes the indicator function on the set 
A x B (see [TJ Theorem 2.6.4, p. 105]). This motivates the following: 

Definition 2.3. Consider a random vector (Xi, . . . , X n , X n+ i, . . . X n+m ) of non- 
negative random variables. Consider N n G B(R") such that P((Xi, . . . , X n ) G 
N n ) = 1. Lei n be a probability kernel from N n to R™ verifying for each Borel sets 
ACN n , BC R™ 

/ dF n+m (x,y)= / dF n (x)ff(B). (3) 

./{xeA.yeB} 7a 

T/ie family {/i x (-), x G iV n } will be said a N n regular conditional distribution 
of (X n+ i, . . . , X n+m ) given (Xi, . . . , X n ). A family of m-dimensional random vec- 
tors {(Z* +1 , . . . , Z* +m ), x G N n }, such that for each x G N n , the distribution 
of (Z* +1 , . . . , Z* +m ) is given by yU x (-); will be said a N n - distributional version of 
(X n+1 , . . . , X n+m ) given {X u . . . , X n ) . 

Remark 2.4. The following remarks are in order: 

1. If in the previous definition the distribution of {X\, . . . , X n+m ) is discrete 
(resp. absolutely continuous), the set N n can be taken as the set of points 
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x G R™ such that (Xi,...,X n ) = x has strictly positive probability (resp. 
strictly positive density function), and the conditional distributions, defined in 
the usual way, satisfy (Ej). Another important example is when (X x , . . . ,X n ) 
and (X n+ i, . . . , X n+m ) are independent. In this case we can take N n = R™ 
and /i x (- ), constant for all x ; being the distribution of (X n+1 , . . . , X n+m ). 

2. The random vectors (Z* +l , . . . , Z* +m ) in Definition \2.3\ are introduced for no- 
tational convenience (actually, they don't need to be defined on the same prob- 
ability space, for different values of x). In fact, the properties we need to use 
(stochastic monotonicity, for instance) will apply to the vector {Z* +l , . . . , Z* +m ), 
with x fixed and have to do with the probability distribution of this vector, not 
with the paths of the process. 

3 DFR property under decreasing conditions of 
the inter arrival times 

In this Section we present our first result concerning the DFR preservation property 
in a renewal process, making use of the conditional distributions of the inter renewal 
epochs. Let {N(t) : t > 0} be a general counting process with inter renewal epochs 
(X„) n=0) i i .... Let T be a DFR random time independent from the process, whose 
survival function is denoted by Ft- We know that if T and X\ don't have as a 
common discontinuity point, then (cf. [HI Lemma 3.2]) 



P(N(T) > n) = E(F T {S n )) 



n = l,2,... 



(4) 



and inequality ([!]) can be rewritten, for each n as: 



E^FriX,)) < E[F T (S 2 )} 



n = 0, 



(5) 



E 2 (F T (S n+1 )) < E(F T (S n ))E(F T (S n+2 )), n = l,2,.... 



(6) 
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In the next result we give sufficient conditions in order to verify ([6]) when we can 
find distributional versions of (X n+ i, X n+2 ) given (Xi, . . . ,X n ) satisfying (jSJ). 

Proposition 3.1. Let {N(t) : t > 0} be a general counting process with inter 
renewal epochs (X n )n=o,i,...- Let T be a DFR random time independent from the 
process. Let n be a fixed natural number and assume that there exists a N n distribu- 
tional version of of {X n+1 , X n+2 ) given (X x , X n ) (say {(Z* +1 , Z* +2 ), x G N n }) 
satisfying, for every DFR survival function H with H(0) = 1, that 

E 2 [H(Z* +1 )] < E\H{Zl +1 + Z* +2 )], x G N n . (7) 

Then, 

E 2 (F T (S n+1 )) < E(F T (S n ))E(F T (S n+2 )). (8) 



Proof: In order to check inequality ([H]), consider x := ( ) G N n and call 

s n := X\ + . . . + x n . As F T is a DFR survival function, then the survival function 
H defined as 

77 W = %±£>, ,>o (9) 

F T (s n ) 

verifies the DFR property (see [TS], p. 118])) and, as H(0) = 1, condition (JTj) holds 
true. If we call /i x to the common distribution of (Z* +1 , Z* +2 ), we can write ([7]) as 

1/2 



H(x n+1 )d^(x n+u x n+2 ) < / H(x n+1 + x n+2 )d/j x (x n+1 ,x n+2 ) 
which, recalling ([9]), turns out to be 

/ F T (s n + x n+1 )d/j*(x n+1 , x n+2 ) 

1/2 

< Ft 2 (s n ) ( / ^ F T (s n + x n+ i + x n+2 )dfi x (x n+1 , x n+2 ) ) . ( 10) 
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flr(x) := / F T (s n + x n+ i + x n+2 )dfi*( ), xGiV n . (11) 

</r2_ 

Integrating (TIP]) with respect to F n , the distribution function of (Xl, . . . ,X n ), we 
have 

/ dF n (x) f F T (s n + x n+l )df(x n+l ,x n+2 ) < [ dF n (x)F^ /2 (s n y /2 (x). (12) 



Applying Cauchy-Schwartz inequality to the term in the right-hand side of (112"]) we 
obtain 

/ dF n (x)F^( Sn y /2 (x) <( [ dF n (x)F T (s n )) ' ( [ dF n (x)g(x.)) ' . (13) 

JN n \JN n J \JN n J 

Thus, by ( !T2|) . (|T3"|) and recalling (TTTj) we can write 

< / dF n (x.)F T (s n ) / dF n (-x) / F T (s n + x n+i + x n+2 )dfi x ( 

The previous inequality shows (jHD □ 

Two results concerning the preservation of the DFR property were given in [6]. 
The first result (See Section 3 in [6]) involves the association (see Definition 15.21 in 
Section 5 below) and stochastic decrease of the interarrival times. In the second 
result (see Section 5 in [6J) a different set of conditions were given. These conditions 
also involve a decreasing behaviour of the interarrival times (together with some 
other technical conditions). We now introduce a different way to check (jSJ) and (El), 
assuming a stronger decreasing condition on the interarrival times than in the two 
previous results, but with no more additional conditions. This property will allow us 
to check the d-DFR condition for generalized renewal processes under more general 
assumptions than in [HI Section 5]. 
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Proposition 3.2. Let (Xi,X 2 ) be a nonnegative random vector. Assume that Ft is 
a DFR survival function. Assume also, that there exists a N\ distributional version 
of X2 given X% (say {Z 2 X , x\ G N{\), such that Z 2 l <st X\, for all x\ G N\. Then, 

E 2 [F r (X!)] < E[F T (X 1 + X 2 )]F T (0). 



Proof: First of all, if Ft is is a DFR survival function, then Ft{- )/i*V(0) is a 
DFR survival function and, therefore, NWU (see [THJ p. 181]). Thus, 

M±*i > Mkw, M >„. 

F T (0) - (Ft(0)) 2 
Then, choosing z = X\ and t = X2 and taking expectations 

F T (0) J B(F T (X 1 + X 2 )) > E{F T {X 1 )F T {X 2 )) (14) 

Now, consider the family {Z^ 1 , x\ G iVi}, and call {fJ-i 1 , X\ G iVi} to its associate 
family of distributions. Recalling (j2J), we can write 

E{F T {X X )F T {X 2 )) = [ F T (x 1 )dF 1 (x 1 ) [ F T (x 2 )d^(x 2 ) 

J JVi JR+ 

= [ F T (x 1 )E(F T (Z 2 a;i ))rfF 1 (x 1 ). (15) 
As we have in this case that Z^ 1 <st X\, then 

E(F T (Z^)) > S(F T (Xi)). (16) 
Thus, we deduce from (JT5]) and ffTB"]) 

E(F T (X!)F T (X 2 )) > ^(FtM). 

This inequality together with ( IT41 prove Proposition 13.21 □ 

Our aim now is to show the preservation of the DFR property making use of 
Propositions 13.11 and 13 . 21 This involves to find a distributional version of (X n+ %, X n+2 ) 
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given (Xj_, . . . , X n ) satisfying the conditions in Proposition 13.21 However, in many 
situations, unidimensional conditional distributions are given in a natural way, that 
is, we know the distribution of X n+ i given (Xi, . . . , X n ), n = 1, 2, . . . . In the dis- 
crete, or absolutely continuous case, the construction of bidimensional conditionals 
provided sequentially unidimensional conditionals is well-known. In order to extend 
this formula, recall that if we have a probability kernel fi n +i from N n G B{R r \_) to 
R + , and we have a probability kernel /i n +2 from N n x R + to R + , we can define the 
composition kernel v := fi n+ i (g) fi n+ 2 from N n to R+ as (cf. [TSJ p. 20] 

i%(B 1 xB 2 )= [ dt£(x n+1 ) [ dfi ( ^ +1 \x n+2 ), xGiV n , (17) 

JB\ J B 2 

for any Borel sets B\ G £>(R + ) and B 2 G B(R + ). In the next Proposition we show 
that composition kernels are the natural way to construct bidimensional conditional 
probability kernels. Its proof, although using standard techniques, is rather long, so 
that it is postponed to the appendix. 

Proposition 3.3. Let (X n ) n=0i i v .. be a family of inter renewal epochs. For n fixed, 
call X = (Xi, . . . ,X n ) and assume that {/i^, x G iV n } is a N n regular conditional 
distribution of X n+ \ given X. Assume also that {fi^+i n+1 \ ( x > x n+i) £ X n x R+} is 
a N n+ \ := N n x R + regular conditional distribution of X n+ 2 given (X, X n+ i). Then 
we have: 

a) A N n regular conditional distribution of (X n+1 ,X n+2 ) given X is given by the 

composition kernel v := fi n+ i ® /i n +2 defined in (77] ). 

b) Take a family of bidimensional random variables {(Z* +1 , Z* +2 ), x G N n }, such 

that for each x G N n , the distribution of (Z* +1 , Z* +2 ) is given by z/* ; as in 
part a). For each x G N n , the family x n+i G R+} is a R + regular 

conditional distribution of Z* +2 given Z* +1 . 

Now we are in a position to prove one of the main results of this section. 
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Proposition 3.4. Let {N(t) : t > 0} be a general counting process with inter 
renewal epochs (X n ) n =o,i,...- Let T be a random time independent from the process. 
Assume that: 

a) T is a DFR random variable. 

b) T and X x don't have simultaneously positive mass at 0. 

c) For each n — 1,2,3, ... , there exists a N n distributional version of X n+ i given 

(X lt X n ) (say {Z* +1 , x G N n }) verifying that 

cl) Z? < ST X U Xl eNx, 

c.2) For n = 1,2,..., N n+l CiV„xR + and 

Z%.l n+l) <st Zf l+l , for all (x, x n+1 ) G N n+1 . 

Then N(T) is discrete DFR. 

Proof: The result will follow by checking ([5]) and ([6]). To the first inequality, 
note that assumption cl), allow us to apply Proposition 13.21 thus having 

E 2 (F T {X 1 )) < E{F T {X X + X 2 ))F T (0) < E{F T (X X + X 2 )), 

and therefore fl5]) holds true. As for ([6]), thanks to c.2), we have, for fixed n, the 
family of random variables {Z^ 2 n+1 \ ( x ; x n+x) £ N n+ \}, being a iV n+1 distributional 
version of X n+2 given (X±, . . . , X n +i). We extend it to a N n x R + version by taking 
Z n +2™ +1 = f° r eacri ( x J x n+i) £ N n x R + such that ( x n +i ) G" N n+ i. Then, we 
have 

^i+2 n+l) <st Z* +1 for all (x, x n+1 ) GiV„x R + . (18) 

This version, together with {Z% +1 , x G N n }, the N n distributional version of X n+ i 
given (Xi, . . . , X n ), provides us, applying Proposition 13.31 a), a N n regular condi- 
tional distribution of (X n+ i, X n+2 ) given X, with its corresponding distributional 
version {(Z% +1 , Z* +2 ), x G N n }. 
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By Proposition 13.31 b) {Z^ n+ , x n +i € R + } is a R + version of Z* +2 given 
Z*+i, which verifies (fl8l) . so that we can apply Proposition 13.21 to the random vector 
(Xi,Xz) ='■ (Z* +1 , Z* +2 ) and any arbitrary DFR survival function H such that 
if (0) = 1, thus having, 



Thus, conditions in Proposition 13.11 are verified, and (pj) holds true. This completes 
the proof of Proposition 13.41 □ 

4 Applications to maintenance models 

Maintenace models describe the behaviour of a mechanism subject to repair. When 
the mechanism fails a repair is performed. This repair can restore the system to 
its initial state (perfect maintenace) or to the condition as it was before the failure 
occurred (minimal maintenace). Also, it can leave the system in an intermediate 
state between both of them (imperfect maintenace), or even in an operational state 
worse than just before the failure (worse maintenace). This kind of models were 
introduced by Kijima and Sumita (see p2]). F° r a recent review of these models we 
refer to J8]. Among them, the so-called virtual age models, associate, to each time 
t, the virtual age of the system v(t), which represents the operational stage of the 
system after successive failures and its corresponding repairs, in such a way that, if 
G represents the lifetime of a new unit, then the lifetime of a unit with virtual age 
v is given by 



In virtual age models the underlying renewal process is defined through the inter- 
arrival times (A n ) n=0j i i ... representing the time elapsed between consecutive failures 
of the mechanism. We also have the sequence (Ki)n=o,i,... representing the virtual 
age of the system just after the n-th repair (we assume there is no delay due to the 




(19) 




v > 0. 



(20) 
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repair). The virtual age of the system usually depends on (X n ) n=0i i i ... and a fam- 
ily of nonnegative random variables (A n ) n= o ) i,..., representing the degree of repair 
performed at n. The value A n = represents the best possible repair of the model 
under consideration, whereas A n — 1 stands for minimal maintenace (the repair 
leaves the system as it was just before the failure). For instance, Kijima (cf. [16]) 
considered two models, in which the virtual age is defined as follows: 

V n = V n -i + A n - X n , Kijima type I, (21) 
V n = A n - (K-i + X n ), Kijima type II. (22) 

In general, the virtual age is defined, for n — 1, 2, . . ., as 

V n :=*(V n - 1 ,X n ,A n ) = * n (Xi,...,X n ,A 1 ,...,A n ), V = 0, (23) 

in which \1/ is a suitable function. Once the virtual age at instant n is known, the 
survival function of the next repair is evaluated as follows. 

P{X n+l > z\V n = v) = G(z\v), v>0. (24) 

In the particular case of (A n ) n= o,i,... being deterministic nonnegative numbers, it is 
clear by ( 123"]) and fT24"|) that a i?" distributional version of X n+ i given Xi, . . . , X n is 
given by a family of random variables x G -R™}, whose reliability function 

is given by 

P(Z* +1 >z) = G(z|tt„(x, A x , . . . , A n )), z>0. (25) 

Of particular interest is the case in which A n — q, n = 1, 2, . . ., for a given < q < 1. 
Then, we obtain the generalized renewal processes considered by Kijima and Sumita 
in (HZ]). 

In the following Proposition we prove the d-DFR property for N(T), for general 
Kijima type I models. The sequence (^4 n )n=o,i,... can be random, although assume 
that, for each n, A n+ i is independent of X\, . . . ,X n . This assumption is also made 
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in the original paper of Kijima [16], along with the assumption that the sequence 
(A n ) n= o,i,... is formed of independent random variables, as well as < A n < 1. We 
don't need to make these last assumptions, as the method of proof is practically the 
same. This result can be useful when an age replacement is used up to a random 
time T, in which a preventive maintenance can be performed. 

Proposition 4.1. Let {N(t) : t > 0} be the counting process associated to a Ki- 
jima type I repair model as in l[21\) . Assume that the sequence of random repairs 
(A n ) n=0) i i ... satisfies that A n+ i is independent of (Xi, . . . , X n ), n = 1, 2, . . .. Let T 
be a random time independent from the process. If the following conditions hold 
true: 

a) T is a DFR random variable, 

b) Xi is increasing failure rate. 

Then, N(T) is discrete DFR. 

Proof: First of all, we will consider that the values in sequence (^4 n )n=o,i, -- * s 
formed of deterministic numbers. The result will follow applying Proposition 13.41 
As Xi is IFR, then it doesn't have positive mass at (see ( [TH|, p. 109])), so that 
condition b) is satisfied. As for c), note that, taking into account (I2T1) . in this case 
ipn m (|23l) has the form 

* n (Xi, ...,X n ,A u ...,A n ) = AiXx + . . . + A n X n . (26) 

As Xi is IFR, then G(- \y), as in (|20|) is decreasing in y, and therefore, taking into 
account ff24l) and (|26l) . we have 

P(Z? >z) = G(z\A lXl ) < G(z\0) = P(X 1 > z), (27) 

so that Zf <st Xi, and c.l) is verified, c.2) is checked in a similar way. In this case 
we have 

P(Z^ +1 >z) = G{z\A lXl + ... + A n+1 x n+1 )} 
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< G(z\A lXl + ... + A n x n )\ = P(Z* +1 > z) 

and Z*+2 +1 <5T Z£+X' Thus, conditions in Proposition 13.41 are verified when the 
sequence (Ai)n=o,i,... is formed of a sequence of constant values and in this case, 
N(T) is d-DFR. Now, assume that the sequence (Ai)n=o,i,... is formed of random 
variables. Fix n a natural number. In order to show (CQ) we will condition to the 
random vector n := (Ai, . . . , A n+ i), that is 



P(N(T) >n + i)= / P(N(T) >n + i\@ n = 9)dF &n (9), i = 0, 1, 2. (28) 

For fixed 9 := (ai, . . . , a n+ i), consider a Kijima type I model {N e (t) : t > 0} 
as in ( 12T]) . in which (A n ) n= o,i,... verifies = aj, i = 1,2, ...,n + 1 (the rest of 
values in the sequence can be defined arbitrarily). Note that, taking into account 
(H]), P(N(T) > n) depends only on (Xi, . . . ,X n ), whose distribution is specified in 
terms of X\ and the virtual ages V\ = A\Xi, . . . , V n -\ = A\Xi + . . . + v4 n _iX n _i), 
so that using the fact that each Ai is independent of the process of past lifetimes 
X\, . . . , X n _i, 

P(N(T) >n + i\B n = 6) = P(N e (T) =n + i), 2 = 0,1,2, (29) 
and fl2HD and ([29]) gives us 

P(N(T) > n + 1) = f P(N e (T) >n + l)dF &n (9), (30) 

JR™ +1 

By the previous equality and taking into account that {Ng(t) : t > 0} is a Ki- 
jima type I model with deterministic repairs (and therefore, verifying the d-DFR 
property) we obtain 

P(N(T) >n + l)< [ P*(N e (T) >n + 2)P*(N e (T) > n)dF &n (9) 



'R" 

l 
2 



< / P(N e (T) >n + 2)dF &n {9)\ / P(Ng(T) > n)dF &n (9) 



Rn+1 / \ /T>n + 1 
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Now, we use the previous inequality together with ([25]) and (|29|) to conclude 

P{N{T) > n + 1) < P5 (N(T) >n + 2)P^ (N(T) > n) □ 

Remark 4.2. T/ie previous result is more general than the one given in Corollary 
5.1.], in two senses. Firstly, we allow the repairing sequence to be random and 
secondly, the initial repair Xi is more general. The question is that in we used 
an analogous proof than the one considered in J&j/ to obtain IFR preservation results, 
and for this property stronger conditions on the arrival epoch are needed. 

Also, it is interesting to point out that, under the same conditions as in Propo- 
sition\4^J\ Kijima type II models don't verify the DFR property for N(T). To show 



that, consider a Kijima type II model, in which X± is increasing failure rate, and 
A\ — \, A2 — (the problem here is that the second interarrival time decreases with 
respect the first one, but, as the second repair is perfect, the process at this point 
'restarts as new'). Consider that T is an exponential random variable of mean 1 
(that is, Ft{x) = e~ x , x > 0). Firstly, we observe in l[21\) . |HP that the virtual age 
V\ in Kijima type I and Kijima type II coincide (as Vq = 0) and, taking into account 
[2l\ ) and Proposition \3.2[ the discrete DFR property for n = 1 is satisfied, that is, 

E 2 [e- Xl ] < E[e- {Xl+x ^]. (31) 

On the other hand, as A 2 = 0, we see from ( [H| ) that V2 = 0, so that X 3 =j X\, and 
is independent of Xi and X 2 . Thus, we can write 

E[e- S *} = E[e~^ +x ^]E[e- Xs ] = E[e^ +x ^}E[e- Xi }. 

Assume that the discrete DFR property is true for n = 2. Taking into account 
and the previous expression we have 

E 2 [e - { x 1+ x 2)] < E[e- s *\E[e- x ^ = E[e-^ +x ^]E 2 [ e -^]. 
We conclude therefore that, 

E[ e ~^+ X ^] < E 2 [e~ Xl ]. (32) 
16 



Thus, taking into account ( Tff7]j and (CUP, the d-DFR for n = 1 and n = 2 would 
imply 

E[e -(X 1+ X 2)] = E 2 [e - Xl ^ 

which is not, in general, true. Take for instance X% to be uniformly distributed on 
(0, 1). Then, we can check easily that 

P(Z? > z) = G(z\ Xl ) = 1 ~ Xl ~ Z , 0<z<l- Xl 

1 — X\ 

and thus, the distribution of is uniform on (0, 1 — x\). Then, we can write 

rl rl-xi -X2 rl rl 

E[e- {Xl+X2) ] = / d Xl e~ Xl / dx 2 = / d Xl e~ Xl / due~ {1 - Xl)u , 

Jo Jo 1-xx J Jo 

which is clearly greater than 



E 2 [e~ Xl ] = [ d Xl e~ Xl [ due~ u . 
Jo Jo 



5 The DFR property under association of the in- 
terarrival times 

In [B], Proposition 3.1 (a) sufficient conditions for two random variables in order to 
satisfy jSJ) were given. These conditions involve the concept of association (cf [TT]). 
We now recall this concept. 

Definition 5.1. A random vector X := (X±, . . . ,X n ) is said to be associated if 

Cov(f(X),g(X)) = E[f(X)- g(X)\ - E[f(X)]- E[g(X)} > (33) 

for all increasing functions f, g : R n — > R such that the expectation exists. 

On the other hand, we know from [6J that 

Proposition 5.2 ( [6J, Proposition 3.1(a)). Let (X n ) n= o,i t ... be an arbitrary sequence 
of interarrival times and T a random time independent of them. Assume that T is 
a DFR random variable. If (X\,X2) is associated and X 2 <st X\, then 

E 2 \F T {X X )\ < E[F T (S 2 )]F T (0). 
17 



As an immediate consequence of the previous result and Proposition 13. II we have 
the following 

Corollary 5.3. Let {N(t) : t > 0} be a general counting process with inter renewal 
epochs (X n ) n= o i t .„. Let T be a random time independent from the process. Assume 
that 

a) T is a DFR random variable 

b) T and Xi don't have as a common discontinuity point. 

c) (Xi,X 2 ) is associated and X 2 <st X\. 

d) For fixed n there exists an N n distributional version of (X n+ i,X n+2 ) given 
(Xi, ...,X n ) (say {(Z* +1 , Z* +2 ), x G N n } satisfying: 

(Z* +1 , Z* +2 ) is associated and Z* +2 <st Z* +1 , for all x G N n . 

Then, N(T) is discrete DFR. 

Proof: We will show that, under these hypothesis, §5$) and are verified. In 
fact, Condition c) and Proposition 15.21 imply immediately (j^D- To prove Q, fan 
and consider H a DFR survival function with H(0) = 1. Condition d), together 
with Proposition 15.21 applied to (Xi,X 2 ) := Z* +2 ) implies that 

E 2 [H(Z: +1 )} < E[H{Zl +l + Z* +2 )\. 

This, together with Proposition 13.11 proves (fU), and therefore the d-DFR property 
for N(T). □ 

Remark 5.4. In f^j, Proposition 3.2. the same condition (c) was required, but it was 
also required the association of the random vector (S n , X n+1 , X n+2 ), for n = 1, 2, . . .. 
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It is not hard to find a renewal process satisfying the conditions in Corollary \5.3\ 
but not verifying the previous association condition. For instance, take Y a random 
variable with support on [1, oo), and let (W n ) n= o,i,... be a sequence of independent and 
identically distributed non-negative random variables independent ofY and verifying 
Y >gji W\. Define the inter-renewal epochs as follows 



X X = Y X 2 



W\ and X n = — — — , n 



3,4, 



(34) 



It is easy to check that conditions in Corollary \5.3\ are verified. However, (X x , X 2 ,X 3 ) 
is not associated. To this end, take Y verifying that 



E[Y]-E 



1 

Y 



> 1 



(35) 



(this inequality is satisfied, for instance, ifY = l + Ber(p) with Ber(p) a Bernoulli 
random variable with probability of success < p < 1). Now, consider that Y and 
Wi have finite mean. We will show that (X\,X 3 ) is not associated. To this end, 
take the increasing functions 

f(x, y) = x, g(x, y) = y, (x, y) G R 2 . 



Taking into account (3$) , ([25P, the fact that Y and W 2 are independent and that 
E[W 2 ] > 0, we have 



E[f{X x ,X 3 )g{X x ,X z )\ 



E[W 2 ] < E[Y]E[W 2 }E 
~W 2 



1 

Y 



E[Y]E 



Y 



E[f{X 1: X z )].E[g{X u X 3 )\ 



Thus (X X ,X 3 ) is not associated and therefore (X x , X 2 , X 3 ) is not associated (see 
fTPj p. 123]). In general, the association of random variables does not imply the 
conditional association and vice versa (see [23] for some examples)) so that the 
conditions given in Corollary 15.31 and in [6], Proposition 3.2. are different. The 
practical advantage in this last result is that we don't need to build the conditional 
distributions of the interarrival times. 
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A Appendix 

In this Appendix we include the proof of Proposition 13.31 

Proof of Proposition 13.31 Let {z/*, x G N n } be the composition kernel, as 
defined in part (a). To show that this kernel is a N n regular conditional distribution 
of (X n+1 ,X n+2 ) given (Xi, . . . , X n ), we need to show for each A G B(N n ) and 
B G B(R 2 + ) (recall ©) 

/ dF n+2 (yi, X n+ i, X n+ 2) = 

)• (36) 

JAxB J A JB 

Firstly, we will verify the previous equality for Borel sets A G N n , and B of the form 

B :=B 1 x B 2 , Bi G B(R+), i = 1, 2 (37) 

and by a classical extension argument we will conclude that (!36|) is also veri- 
fied for Borel sets B G R+. To show (l36l) . note firstly that, by assumption, 



+i n+1 5 (Xj^n+i) & N n x R + } correspond to a iV n x R + regular conditional 
distribution of X n+2 given (X, X n+1 ), so that we can apply ([3]) to this version, and 
write for B as in fl37|) 

/ dF n+2 (x,rc n+ i,2„ +2 ) 

J AxBix B2 



JAxB! Jb 2 



(3* 



Consider now the function 



= / rf/i^ n+l) (x n+2 )- 1a(x)- l Bl (x n+1 ), 
Jb 2 

As by hypothesis {/x*, x G iV n } correspond to a iV n distributional version of X n+1 
give X, we can apply ([2]) to the previous function and write the right-hand side term 
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in (BHD as 



/ dF n+ i(x,x n+ i) / d^l x { l+l) (x n+2 ) = Eg(X,X n+ i) 

J AxBi JB 2 

= / dF n (x) / d/J,l(x n+1 )g(-x,x n+1 ) 

= / dF n (x) / d//*(x„+i) / d^+{ l+l) (x n+2 ) 

J A J B\ JBo 



), (39) 

A JB 1 xB 2 

where, in the last equality we have used (fT7|) . Therefore, ( 13"S]) and (j3U|) show f[3"6T) for 
a fixed Borel set A G and sets of the form (13"T|) . Using a classical monotone class 
argument (cf. [121 Thm. 1.1., p. 2]) (|36|) holds true for all Borel sets B G R+. This 
means that the family {v*, x G A^ n } satisfies Definition 12.31 thus being a iV n regular 
conditional distribution of (X n+ i, X n+2 ) given X, and part (a) is completed. Part (b) 
follows easily by (a). Let {(Z* +1 , Z* +2 ), x G A^ n } be the family of random variables 
as given in (b). First of all, note that for fixed x G N n , {/4j+i n+1 \ %n+i £ is a 
probability kernel from R + to R + . Secondly note that ( TT7j) gives (J3J) for (Z* +1 , Z* +2 ) 
in terms of the marginal of the first variable (Z% +1 ) and {/^h^l™ > x n +i G R+}, so 
that Definition 5 is verified and we conclude that the family {/4 l +i n+1 , £n+i G R+} 
is a R + regular conditional distribution of Z* +2 given Z* +1 , as claimed. 
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